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PREFACE 


This  meeting,  sponsored  by  the  Joint  Army-Navy-NASA-Air  Force 
(JANNAF)  Performance  Standardization  Subcommittee,  is  held  annually 
to  promote  the  exchange  of  technical  information  among  governmental, 
industrial,  and  academic  scientists  concerned  with  the  experimental 
measurement,  analytical  prediction,  correlation,  extrapolation, 
and  flight  confirmation  of  the  performance  of  liquid  and  solid 
propulsion  systems. 

The  meeting  was  held  February  14-15,  1980,  at  the  Red  Lion  Inn, 
in  Sacramento,  California. 


Henry  F.  Hege 
CPIA  Representative 
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MEETING  SUMMARY 


I.  The  thirteenth  meeting  of  the  Performance  Standardization  Subcommittee 
was  held  14-15  February  1980,  at  the  Red  Lion  Inn,  Sacramento,  CA.  Welcoming 
remarks  were  made  by  Dr  Daweel  George,  AFRPL,  Subcommittee  Chairman.  Mr  M.  Ditore 
ASPC,  served  as  Program  Chairman  and  coordinated  meeting  arrangements.  A  list  of 
attendees  and  the  meeting  agenda  are  presented  in  Appendices  1  and  2. 

II.  PRESENTATIONS.  Reproductions  of  the  slides  and  handout  material  from  the 
meeting  are  given  as  Appendices  1  through  25  in  the  order  that  they  appeared 

on  the  meeting  agenda. 

III.  GENERAL  MEETING  MINUTES. 

A.  The  Solid  Performance  Program.  The  current  version  of  the  program, 
known  as  the  interim  version,  is  operational,  and  documentation  is  forthcoming. 
Program  objectives  of  predicting  solid  rocket  motor  delivered  specific  impulse 
within  +-0.5%  and  thrust  and  total  impulse  with  +-3%  are  being  met  for  the  test 
cases;  however,  a  comprehensive  program  validation  and  verification  are  yet  to 
be  accomplished.  The  Interim  Version  SPP 

(1)  can  handle  different  size  particles  with  particle  impingement, 

(2)  includes  a  new  particle  size  model  and  drag  law,  a  new  combustion 
efficiency  model,  and  a  new  throat  erosion  model; 

(3)  provides  extensive  plots  and  an  expanded  summary  page 

(4)  Fully  Coupled  Transonic  analysis  (FCT). 

Additional  work  is  planned  for  the  following  modules: 

(1)  Grain  Design  and  Ballistics 

(2)  One  Dimensional  Kinetics 

(3)  Two  Dimensional  Two-Phase  Flow 

(4)  Turbulent  Boundary  Layer 

(G.  Nickerson) 

A  status  report  on  the  FCT  subprogram  indicates: 

-  that  for  REGION  I,  (throat  entrance  region)  the  calculations  work 
well  and  are  reasonably  reliable; 

-  that  for  REGION  II,  (around  backside  of  nozzle  nose  cap)  the 
solutions  are  expensive  and  difficult  to  obtain. 

These  results  indicate  that  the  approximate  transonic  analysis  should  be  used 
for  parametric  studies  and  the  FCT  for  final  performance  prediction.  (D.  Coats) 

A  variety  of  Kinetic  Rate  data  was  presented  for  review.  (D.  Coats, 

M.  Salita) 

The  SPP  approach  to  combustion  efficiency,  involving  dimensionless 
variables,  was  presented.  The  erosive  burning  options  of  the  SPP  were  presented, 
along  with  the  various  modelling  efforts  in  progress;  and  the  problem  of  scale-up 
was  discussed.  (R.  Hermsen) 

The  Grain  Design/Internal  Ballistics  Module  (GDB)  contains  four  sub¬ 
programs:  three-dimensional,  two-dimensional  and  axisymmetric  grain  configurations 
and  ballistics.  The  various  features  of  the  GDB  and  some  examples  of  inputting 
various  motor  grain  geometries  were  presented.  (Appendix  9  ).  (J.  Lamberty) 

The  latest  version  of  the  SPP,  the  "Interim"  SPP,  is  available  now,  upon 
written  request  to  D.  George,  AFRPL.  Requestors  of  the  program  should  enclose  a 
computer  tape  and  state  necessary  tape  specifications.  The  requestor  will  be 


provided  with:  a  copy  of  the  program,  the  thermochemical  tape,  sample  case 
input  and  output,  and  a  manual.  Users  are  requested  to  offer  feedback  to 

D.  George.  The  next,  "1980,"  version  of  the  SPP  is  anticipated  to  be  available 
in  the  late  fall  1980.  (D.  George) 

B. 

A  descripti 

models  was  presented.  Attention  was  focused  on  mass  transfer.  Presently,  this 
mechanism  is  based  on  a  single  component  species  (A102H)  and  should  be  based  on 
multispecies.  At  present,  there  is  no  adequate  set  of  screened  reactions  to  go 
with  mass  transfer.  It  was  recommended  that  additional  studies  should  be  con¬ 
ducted  before  including  0D3P  in  the  SPP  methodology.  (K.  Hunter) 

C.  Experimental  particle  impact  data,  based  on  experiments  now  being 
conducted  for  the  first  time,  were  presented  for  review.  Using  a  droplet 

enerator,  it  is  feasible  to  obtain  experimental  impact  data  in  nozzle  flows. 

Z.  Chiba) 

D.  The  Space  Motor  Combustion  Spin  Effects  include  burnrate  increase, 
increased  slag,  and  change  in  Isp.  Motor  analysis  uncertainties  include 
particle  density,  particle  size  distribution,  temperature  gradient  across  the 
particle  on  the  surface,  effect  of  gas  flow  and  temperature  after  particle 
deposition  on  inert  surface,  and  particle  propulsivity.  (W.  Brundige) 

E.  The  effect  of  expended  inerts/slag  on  performance  may  be  taken  into 
consideration  by  utilizing  the  three  proposed  definitions  of  Isp,  expressed 

as:  the  propellant  Isp,  the  delivered  Isp,  and  the  effective  Isp.  (j.  Lamberty) 

F.  A  description  of  the  Burning  Rate  Anomaly  Factor  thought  to  be  due 
to  propellant  flow  during  casting  was  presented.  Common  terminology  includes 
BARF,  hump  effect  and  mound  effect.  (T.  Kirshner) 

G.  The  Particle  Impingement  Modelling  discussion  involved  a  review  of 
work  from  several  years  ago,  and  included  a  call  for  more  current  work  in  this 
area.  The  smaller  particles  coated,  and  the  larger  particles  eroded  the  nozzle. 
Impingement  has  an  important  effect  on  thrust  losses  and  must  be  considered 
when  designing  an  optimum  contour.  (W.  Daines) 

H.  The  Nozzleless  Performance  Program,  currently  under  development  at 
ARC,  will  address  the  phenomena  of  erosive  burning,  grain  deflection,  and 
combustion  efficiency.  Nozzleless  motors  operate  at  65-85%  of  theoretical  Isp— 
increased  efficiency  is  obtainable  by  using  propellants  with  low  pressure 
exponent.  In  the  erosive  burninq  model,  scaling  effects  must  be  taken  into 
account.  (M.  Procinsky) 

I.  SUBCOMMITTEE  STATUS.  Annual  report,  accomplishments,  current  task 

areas,  planned  activities  were  presented.  Appendix  19.  (D.  George) 

J.  Application  of  the  constant  fractional  lag  concept  to  kinetics  and 
mass  transfer  mechanisms  was  suggested.  The  premise  is  that  from  the  nozzle 
Inlet  to  the  throat  various  species  lag  by  some  constant  fraction  their  equili¬ 
brium  value  at  the  throat.  This  approach  would  substantially  reduce  computer 
run  time.  The  throat  values  would  provide  the  initial  conditions  for  full 
computation  In  the  supersonic  region  of  the  nozzle  where  mechanism  behavior  has 
a  significant  effect  on  performance.  A  study  should  be  conducted  to  determine 
the  effect  on  performance  prediction  accuracy  by  application  of  this  concept  In 
the  nozzle  Inlet  region.  (S.  Cherry) 
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One-Dimensional,  Three-Phase,  Flow  Reacting  Gas  with  Mass  Transfer, 
on  of  the  various  particle  size  change  mechanisms  which  this  program 


K.  There  is  a  need  to  establish  a  standard  set  of  Thermochemical  and 
Heats  of  Formation  data  and  a  methodology  for  distributing  and  keeping  them 
current.  (D.  George) 

L. .  Efficiency  Definitions  were  offered  as  a  set  of  equations.  It  was 
stressed  that  consistency  of  prooellant  specie  input  within  equations  and 
between  programs  is  important.  The  "multi"  Isp  concept  was  discussed,  i.e., 
delivered  Isp,  propellant  Isp,  and  effective  Isp.  (D.  Coats) 

M.  Boundary  Layer  Analysis.  An  improvement  of  Aerotherm’s  MEIT  code 
was  described,  whicn  is  consistent  with  TBL,  but  has  added  capabilities.  The 
MEIT  shape  factor  should  be  improved.  The  code  can  handle  rough  or  transpired 
nozzles,  which  TBL  cannot.  (M.  Salita) 

Aerojet  is  conducting  a  boundary  layer  analysis  using  a  finite  difference 
method,  as  opposed  to  an  integral  method,  which  was  developed  by  Dwyer, 

University  of  California,  Davis,  and  comparing  results  with  TBL.  (M.  Ditore) 

N.  A  method  of  Vertical  Force  Test  Measurement  was  presented,  with  the 
observation  that  thrust  vs  time  comparisons  differ  from  horizontal  to  vertical 
testing.  (R.  B.  Runyan) 

O.  Various  aspects  of  Digital  Filtering  Techniques  were  described, 
including  examples,  techniques,  and  the  effects  of  filtering.  The  caution  was 
offered  to  beware  of  the  interaction  of  test  equipment  with  an  experiment. 

Weighted  averages  may  be  used  to  smooth  equipment  oscillation  effects.  (R.  Little) 

P.  A  Three  Dimensional  Time  Dependent  Analysis  was  described.  The 
application  involved  a  rotating  nozzle,  rotating  coordinate  systems,  and  a 
chamber  calculation  with  oscillating  pressure  waves.  (J.  Hoffman) 

Q.  A  Static  Test  Panel— C.  Beckman,  M.  Ditore,  E.  Landsbaum— was 

assigned  to  investigate  performance  measurements,  to  examine  what  performance 
parameters  are  currently  being  measured  versus  what  parameters  should  be 
measured.  The  need  was  stressed  for  performance  analysis  and  tests  to  present 
data  in  similar  terms.  See  Action  Item  #5,  p.  5  .  (D.  George) 


R.  Advanced  Performance  Methodology.  The  limitations  of  current  pre¬ 
dictive  capabilities  were  enumerated,  in  part,  and  then  followed  by  open  forum 
recommendations  to  accommodate  anticipated  future  needs.  These  are  the  ability 
to  calculate,  handle  or  more  accurately  predict  the  following: 


-  large  shock  waves 

-  gouging,  uneven  nozzle  erosion 

-  particle  size 

-  particle  impingement  effect 

-  distributed  combustion  effect 

-  more  kinds  of  metals  (Zirconium,  Boron,  Magnesium) 

-  3-D  grains 

-  multi  flow  ballistics 

-  occurrence  of  tin  in  zirconium 

-  high  energy  propellants 

-  temperature  effect 

-  up  to  300:1  area  ratio 

-  nozzle  submergence 

-  Ignition  transient,  thrust  termination 


(J.  Levine) 
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S.  Meeting  Comments.  The  procedure  for  papers  worked  well  this  year, 
and  the  premeeting  material  contained  most  of  the  papers.  During  the  meeting, 
there  should  be  less  In  the  way  of  presentation  and  more  in  the  way  of  discussion. 
The  1981  subcommittee  meeting  tentatively  will  be  arranged  for  by  Chemical 
Systems  Division,  Sunnyvale,  CA,  (R.  Hermsen)  in  February  1981.  (D.  George) 
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ACTION  ITEMS 


PERSONNEL* 


1.  Write  the  procedure  for  handling 
the  analysis  of  expended  Inerts. 

2.  Chair  a  workshop  on  spin  effects 
of  space  motors  (winter  1981). 

3.  Distribute  write-up  on  Standard 
Method  for  Burn  Rate  Determination. 

4.  Establish  a  consistent  data  bank  of 
thermochemical  data  and  heats  of 
formation,  and  provide  a  mechanism 
for  its  distribution. 

5.  Static  Test  Panel— to  examine  what 
performance  parameters  are  currently 
being  measured  versus  what  parameters 
should  be  measured.  Where,  how,  and 
how  many  measurements  are  appropriate? 
How  should  the  data  be  handled  in  its 
reduction  and  interpretation? 

6.  Write  a  procedure  for  handling  slag 
analysis  in  rocket  motors. 

7.  Provide  write-up  on  reporting  and 
format  of  performance  data  from  motor 
test  firings. 


J.  Lamberty 
W.  Brundige 
D.  George 

D.  George,  C.  Selph 

C.  Beckman,  M.  Ditore, 

E.  Landsbaum 

0.  Lamberty,  W.  Brundige 
J.  Lamberty,  W.  Brooks 


*  First  name  on  list  should  initiate  action  and  coordinate  with  others 
listed  to  resolve  the  action  items. 
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13th  PERFORMANCE  STANDARDIZATION  SUBCOMMITTEE  MEETING 


Abbett,  M. 

Acurex-Aerotherm/Mountaln  View 

Beckman,  C. 

AFRPL/Edwards  AFB 

Brooks,  W.  T. 

Hercules /McGregor 

Brundige,  W.  N. 

Thiokol/Elkton 

Chan,  G.  0. 

Aerojet  Strategic/Sacramento 

Chang,  I-Shih 

Aerospace/El  Segundo 

Cherry,  S. 

KVB/Tustin 

Chiba,  Z. 

Acurex/Mountain  View 

Coats,  D. 

SEA/Santa  Ana 
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Aerospace/Los  Angles 

Lasley,  G. 
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Levine,  J.  N. 

AFRPL/Edwards  AFB 

Little,  R.  R.,  Jr. 

ARO/Aroold  AFS 

Lyon,  J.  M. 

NOS/Indian  Head 

Mikkelsen,  C.  D. 

Thiokol/Huntsville 

Mockenhaupt,  J.  D. 

Aerojet  Strategic/Sacramento 

Nickerson,  G.  R. 

SEA/Santa  Ana 

Procinsky,  I.  M. 

ARC/Gainesville 

Radke,  R. 
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Appendix  2:  Meeting  Agenda 


AGENDA 

PERFORMANCE  STANDARDIZATION  SUBCOMMITTEE 
FEBRUARY  14-15,  1980 

RED  LION  INN 
2001  Point  West  Way 
SACRAMENTO,  CALIFORNIA 


Thursday,  February  14 


Del  Paso  Room 


08:30  A.M.  Opening  Remarks 

D.  George,  PSS  Chairman 
AFRPL,  Edwards,  CA 

08:45  A.M.  Improved  Solid  Performance  Program  (SPP) 

G.  Nickerson,  D.  Coats— Software  and  Engineering  Associates, 

Santa  Ana,  CA 

R.  Hermsen,  J.  Lamberty— Chemical  Systems  Division,  Sunnyvale,  CA 

10:00  A.M.  BREAK 

10:30  A.M.  One-Dimensional,  Three-Phase  Flow  Reacting  Gas  with  Mass 
Transfer 

S.  Hunter— KVB  Engineering, Tustin,  CA 

11:00  A.M.  Particle  Impact  Erosion 

Z.  Chiba— Acurex,  Mountain  View,  CA 

11:30  A.M.  Performance  Methodology  Status  Summary 

J.  Hoffman— Purdue  University,  W.  Lafayette,  IN 

12:00  Noon  LUNCH 

01:30  P.M.  Space  Motor  Combustion  Spin  Effects 

W.  Brundige—Thiokol/Elkton,  Elkton,  MD 

02:00  P.M.  Expended  Inerts/Slag  Performance  Procedure 

J.  Lamberty— Chemical  Systems  Inc.,  Sunnyvale,  CA 

02:30  P.M.  BREAK 

03:00  P.M.  Burning  Anomaly  Rate  Factor 

T.  Kirschner— Thiokol /Elkton,  Elkton,  MD 

03:30  P.M.  Particle  Impingement  Modelling/Performance  Effects 
W.  Daines— Hercules  Inc.,  Magna,  UT 

04:00  P.M.  Nozzleless  Motor  Performance  Mechanisms 

M.  Procinsky— Atlantic  Research  Corp.,  Alexandria,  VA 


Friday,  February  15 


Del  Paso  Room 


08:30 

09:00 

09:30 

10:00 

10:30 

11:00 

12:00 

01:30 

02:00 

03:00 

03:30 


A.M.  Annual  Report  and  General  Business 
D.  George,  PSS  Chairman 
AFRPL,  Edwards,  CA 

A.M.  Thermochemical  and  Heats  of  Formation  Data 
J.  Rendleman— AFRPL,  Edwards,  CA 


A.M.  Efficiency  Definitions 

D.  Coats--Software  and  Engineering  Assoc.,  Santa  Ana,  CA 

A.M.  BREAK 


A.M.  Boundary  Layer  Study 

M.  Salita— Thiokol/Wasatch,  Brigham  City,  UT 

A.M.  Vertical  Force  Test  Measurement  for  MX 

Performance  Measurements  and  Data  Reduction 

E.  Turner,  R.  Little— ARO,  Arnold  AFS,  TN 

Noon  LUNCH 

P.M.  Three-Dimensional  Time-Dependent  Analysis 

J.  Hoffman— Purdue  University,  W.  Lafayette,  IN 

P.M.  Advanced  Performance  Methodology 

J.  Levine— AFRPL,  Edwards,  CA 

P.M.  BREAK 

P.M.  Open  Forum,  Items  of  Opportunity,  Future  Plans 

MEETING  CHAIRMAN  AND  ARRANGEMENTS  COORDINATOR: 


M.  Ditore 

Aerojet  Solid  Propulsion  Co. 
Sacramento,  CA 

MEETING  COORDINATOR: 


H.  F.  Hege 

Chemical  Propulsion  Information  Agency 
Laurel ,  MD 


V*>*&rw  • 


IMPROVED  SPP  INTERIM  VERSION  NOW  AVAILABLE  FOR  DISTRIBUTION 


PSS  CONCLUDE  LIQUID  ROCKET  PERFORMANCE  ACTIVITIES 


PERFORMANCE  MEASUREMENTS 
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Appendix  4:  Improved  Solid  Performance  Program 
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SIMPLIFIED  PROCEDURES 
SUMMARY  PAGES 
PLOTTING  CAPABILITY 


MODULAR  STRUCTURE  OF  THE  SPP 


Number  of  subroutines  in  module.  Number  of  subroutines  in  SPP  *  222. 


COMPUTER  PLOTTING 
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(AXIAL  AND  TRANSVERSE) 


Figure  3:  Case  3,  Mach  Number  Contours,  y=1*2,  Source  Flow,  NILPslO 
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Comparison  of  motor  data  with  erosive  burning 
threshold  predicted  by  Beddini.  Region  above 
solid  line  is  expected  to  have  erosive  burning. 
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EFFECT  OF  MOTOR  SCALE  ON  EROSIVE  BURNING 
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Grain  Design  and  Internal  Ballistics  Module 
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LENOIR  ROBILLARD  (DIAMETER) 


iBEPS  APORT  *  1,  ABURN  -  3,  LM  -  10,  R1000  -  .4,  BRMULT  =  1,  $END 


LIFRAM  GRAIN  DESIGN  FILE 


100  OSDB  IB«0,IDESGN*3,ILIST«0,0EnD 

200  IRR  BOOSTER/  4  FIN  FINOCYL  DAK  1/6/80 

300  OCASE  ZR-O. ,0. ,  .08,1.0,  .1,1.3,  .15,1.0,  .2,1.92,  .3,2.2,  .42,2.5, 

400  .52,2.8,  .08,3.08,  .8,3.38,  1.0,3.03,  1.25,3.9,  1.05,4.08, 

300  1.70,4.102,  42.88,4.102,  OEND 

OOO  f INPUT  1  JN- 13,  JN 1*12,  PHAX=45,0END  8  OF  SURFACES  AND  AN6LE  OF  SYMMETRY 
700  0INPUT3  DELTA-28*. 1 ,5*. 01 ,EP3>3. ,IEND  UEB  STEPS  AND  ACCURACY 

800  01NPUT5  NSFC-1 , ITYPE-2, I ANBLE-1 , THETA-O ., PHI-0. , 

900  X«. 79, Y-0., Z*. 5, R=. 79,. 5, UNIN-O., UMAX-45., 

1000  VHIN— 45., VHAX-O., OEND  HEAD  FILLET  TORUS 

1100  0INPUT5  NSFC=2,ITYPE*4,X=0. ,0. ,Y*0. ,0. ,Z«.5,29.8B,R*1 .29, 

1200  UNIN-O. ,UHAX=45.,VHIN=0.,VHAX«29. 38, OEND  MAIN  CYLINDER 
1300  0INPUT5  NSFC-3, ITYPE-2, I ANGLE* I ,X*4.102,Y=0. ,Z=29.88, 

1 400  THETA=90 . , PHI=90 . ,R>3. 1 85, .375 ,UHIN-90 . ,UHAX-1 34 . , 

1500  VHIN-O., UMAX- 180., OEND  FIN  TORUS  1 

UOO  0INPUT5  NSFC-4, ITYPE-2, I ANGLE* 1 ,X— 2.325, Y=0. ,Z*36. 09, THETA-90. , 

1700  PHI»90.,R»5.752,.375,UHIN— 90., UMAX— 44., VHIN-O., 

1800  VHAX-180.,0END  FIN  TORUS  2 

1900  0INPUT3  NSFC-5,ITYPE-4,X-3.425,3.425,Y*0.,0.,Z»34.09,42.88, 

2000  R« , 375 , UNI N-0 . , UNAX-90 . , VHIN-O . , VMAX-8 . , 

2100  OEND  SMALL  FIN  6ROOVE  CYLINDER 

2200  0INPUT5  NSFC=0, ITYPE*3,X*0.,Y». 375, Z«29. 88, IANBLE*1 , THETA-90. , 

2300  PHI-90. .UHIN--1 4. , UMAX-0. ,VMIN— 3.425, UNAX-O., OEND  FIN  SIDE  PLANE 

2400  0INPUT5  NSFC-7, ITYPE-4.X-0. ,0. ,Y-0. ,0. ,Z-29.88,39.93,R-1 .29, 

2500  UHIN-O., UMAX-45. ,VMlN-0., UMAX-11., OEND  SECOND  CYLINDER 

2400  0INPUT5  NSFC-8, ITYPE— 3,X*0. , Y-0. ,Z-39.93, IAN6LE-1 , THETA-O. , 

2700  PHI-0. ,UHIN-1 .29, UMAX- 1 .71 ,VNIN*.912,VHAX®1 .210, 

2800  OEND  NOZZLE  VERT  PLANE  1 

2900  0INPUT5  NSFC-9.ITYPE-2, IANBLE-1 , THETA-O. , PHI-0. ,X*0. , Y-0., Z-39. 93, 

3000  R*1 .71 ,0. , UHIN-O. , UMAX-45. ,VMIN— 90. , VMAX-49 . , 

3100  OEND  NOZZLE  CORNER  TORUS  1 

3200  0INPUT5  NSFC-10,ITYPE-5,X-0. ,0. ,Y«0. ,0. ,Z*39.93,41 .98,R-1 .71 ,3.5, 

3300  UMIN-O. ,UHAX>45. ,OEND  NOZZLE  CONE 

3400  0INPUT5  NSFC-1 1 ,ITYPE«2,IANGLE*I , THETA-O. , PHI-0., X-0. , Y-0. ,Z-41 .98, 
3500  R-3. 5,0., UMIN-O., UHAX-45. ,VMIN— 41., VMAX-O., 

3400  OEND  NOZZLE  CORNER  TORUS  2 

3700  0INPUT5  N8FC-12,ITYPE-4, X-0. , 0. , Y-0., 0. ,Z»41 .98, 42.88, R-3.5, UHIN-O., 
3800  UMAX-45. ,VMIN*0. ,VHAX>2. ,OEND  NOZZLE  CYLINDER 

3900  0INPUT5  NSFC-13, ITYPE— 3, X-0., Y-0., Z-42.8B,IAN0LE-1 , THETA-O., PHI-0., 
4000  UNIN-3. 5, UMAX-4.1 02, VMIN-O. ,VMAX-2. 90, OEND  NOZZLE  END  PLANE  2 

N 


LIFRAM  BALLISTIC  FILE 


100  IGBB  IB=1,IEND 

200  IRRFIN  IBB  BOOSTER  *145  +3SIG  1/4/80 

300  IBALPAH  BT*1 .94,PIN*25. ,PFINAL«1 6. ,WEB*2.81 ,VIN«349.24,ETA*.9348, 

400  R1000«.7938,NEXP*.4,Re*4?.?7,C8TAR«5057.?e,RHOP*.06A1,8AN«1.13,TF*6623., 

SOO  DE*4.94,PA*14.7,PERCNT*10. ,EBTYPE*1 ,NBS*0,NPT*2,NH*1?,BETA*70.,PLT*2,IEND 

600  OPTIHE  XY*0. , .01 ,  .6, .02,  DTAIL*. 005, STIHE-10, BEND 

700  OPLOTTS  IPLOT-1 , 1 ,0,0,PSCALE=.3, JPL0T*1 ,IEND 

800  OEROS  KODE*7,XY*0.,0. ,  .05,0.,  .1,25.,  .15,31.,  .5,32.3,  t.,33.5,  1.5,34.8, 

900  2. 0,36. 1,2. 6, 37. ,2. 95, 35. 3,  3.05,  0.,IEND 

1000  IAPORT  RP0RT*1.29,IEND 

1100  IAPORT  XY*0.,12.52,  2.81,52.81,  IENB 

1200  IBEPS  ABURN*1 ,AP0RT*1,ALFA*. 15,LN*.5,I£ND  HEAD  FILLET  TORUS 

1300  IBEPS  ABURN*2,AP0RT*1, ALFA*. 15, LH*5.88, BEND  MAIN  CYLINDER 

1400  IBEPS  ABURN*2, APORT* 1 , ALFAS.15,LN=5.88,IEND  5  ELEHENT8  OF  1  SURFACE 

1500  IBEPS  ABURN*2,AP0RT*1, ALFA*. 15, LH*5.88, BEND 

1400  IBEPS  ABURN*2,AP0RT*1, ALFA-. 15, LN*5.88, IENB 

1700  IBEPS  ABURN*2,AP0RT*1 ,ALFA*.1S,LH*5.88,!END 

1800  IBEPS  ABURN*3 , APORT *2 , ALf A* . 1 5 , LH=2 . ,IEHB  FIN  TORUS  1 

1900  IBEPS  ABURN*4,AP0RT*2, ALFA*. 15, LH*4., BEND  FIN  TORUS  2 

2000  IBEPS  ABURN*5 , 6 , 7 , APORT  *2 , ALF  A* . 1 5 ,LH* 1 .6, IENB  5  ELENENTS  OF  3  SURFACES 

2100  IBEPS  ABURN*5,6,7,AP0RT*2,ALFA*.15,LR*t.6,!END  FIN  GROOVE  CYLINDER, 

2200  IBEPS  ABURN*5 ,4,7, APORT *2, ALFA*. 15, LN*1. 6, BEND  SIDE  PLANE, 

2300  IBEPS  ABURN*5,4,7,AP0RT*2,ALFA*. 15,LN*1 .6,IEND  AND  TOP  CYLINDER 

2400  IBEPS  ABURN*S,4,7,AP0RT*2,ALFA*.15,LH*I .4, BEND 

2500  IBEPS  ABURN*8,AP0RT*2,LH*.1 ,IEND  NOZZLE  VERT  PLANE  1 

2400  IBEPS  ABURN*9,AP0RT*2,LN*.2,IENB  NOZZLE  CORNER  TORUS  1 

2700  IBEPS  ABURN*10,AP0RT*2,LH*2.05,!END  NOZZLE  CONE 

2800  IBEPS  ABURN*t 1 , AP0RT*2,LH* .2, BEND  NOZZLE  CORNER  TORUS  2 

2900  IBEPS  ABURN*12,AP0RT*2,LH*.9,IEND  NOZZLE  CYLINDER 

3000  HEPS  ABURN*13,AP0RT*2,LN*. 1 , BEND  NOZZLE  PLANE  2 

N 


AX I SYMMETRIC  GRAIN 


Set  up  for  Grain  Design  and  Ballistics  (Sequential  Solution) 


100  BGDB  IB=2,IDES6N*1,0END 

200  TEST  CASE  FOR  AXISYNHETR1C  GRAIN 

300  OCASE  ZR*0,2,  0,5,  1,6,  20,6,  21,5,  21 ,2, SEND 

400  I6RAIN  ZR=0,2,  15,2,  14,4,  15,4,  16,2,  21,2, 

500  ARC=0,0, 1 ,0,0, SEND 

&0o  ISTEP  DELTAZ22* .2, IEND 

700  AXISYHMETRIC  TEST  CASE 

BOO  OBALPAH  DT*1.I3,  PIN=18,  «EB*4,  VIK=251 ,  ETA«.?6,  R1000*.3?5, 

900  NEXP*.25,  RG*56,  CSTAR»5266,RHOP*.063,  6AH-1.I8,  TF*6tOO, 

1000  EIH*B,  BETA'70,  PA-14.7,  PERCNT=10,  EBTYPE*1 ,  HH«3,  NBS*0, 

1100  PFINAl>16,6EMD 

1200  0PT1HE  XY»0,.02,  ,1,.1,  .2, .2,  DTAIL=.02,  DURATN*20,  3TINE«1 ,IEHD 

1300  OEROS  K0DE*7,  XY*0,2,  50,2,0END 

1400  OBRHULT  XY*0,1.5,  4,1.5,0END 

1500  IBEPS  ABURN-1,  APORT-1,  BRHULT*1,  IN-15,IEKD 
1600  OBEPS  ABURN«2,3,4,  APORT*l ,  LH*1 ,!END 
1700  OBEPS  ABURN>5,  AP0RT>T,BRHULT>1 ,OEND 

N 
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700  t GRAIN  TYPE*2,  N*4,Z»57.545tR*2.357fM»1 .382,R1«.21 4,R2*. 195,ETA=60,PSI=30fOEND 
BOO  06RAIN  TYPE«5,R*2. 357, Z»43. 575, XY». 9375,0,  .44281 f .66281 ,  0,.93?5, 

900  DE6«?0,  ARC*-1 ,NE6BRN*0f0f OEND 

1000  06RAIN  TYPE*5,R=2.35r , Z=67. 8, XY=. 9375,0,  .66281 ,.66281 ,  0..9375, 

1100  ARO-1,  NEGIRN*Q,0,0END 


SIDEWINDER  BALLISTICS  FILE 


NOl iVOOl  *3«*  JDVJWtS  ON*  31IHKJH  JIISIVW* 


lo«itrr  Hrpion 


gaB^Bpgas^BSBgaEaghga _ _ 
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aaqmPHMiaHWl 
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Aft  and  Forward  Pressure  Vs  Time 


IGNITER  MODELING 


TIME,  SEC. 
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INTRODUCTION 


GENERAL 

This  is  a  summary  of  the  work  accomplished  on  Air  Force 
Contract  F04611-78-C-0011.  A  computer  program  was  developed 
by  KVB  in  1973  to  analyze  one-dimensional  reacting  gas-particle 
flows  primarily  involving  coal  combustion.  This  program  was 
developed  from  experience  gained  in  prior  work  on  rocket  engines 
and  has  been  re-converted  for  application  to  two  phase  flow 
in  rocket  motors  during  this  contract. 

The  purpose  of  the  current  project  was  to  further  develop 
the  existing  computer  program  to  incorporate  additional  particle 
size  change  mechanisms  that  are  of  importance  in  improving  the 
ability  to  predict  rocket  motor  performance.  Also  the  program 
was  modified  to  simplify  the  input  required  that  is  specific 
to  rocket  motor  performance  prediction. 

All  analytical  work  was  performed  at  KVB,  Inc.,  Tustin, 
CA.  Computer  programming,  coding,  and  computation  was  performed 
by  KVB  with  the  assistance  of  Software  and  Engineering  Associates 
(SEA) ,  Santa  Ana,  CA,  acting  as  consultants  to  KVB.  Dr.  Cye 
Waldman  served  as  a  consultant  to  KVB  performing  the  literature 
survey,  and  selection  and  development  of  several  of  the  size 
change  models. 

OBJECTIVE 

The  objective  of  this  program  is  to  improve  the  ability 
to  predict  the  effect  of  two  and  three  phase  flow  losses  on 
the  performance  of  solid  propellant  motors.  The  improved  predic¬ 
tion  capability  is  obtained  by  the  development  of  the  analytical 
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capability  to  predict  particle  and  droplet  size  change 
mechanisms  in  a  reacting  solid  rocket  nozzle  flow. 


SCOPE 


Five  program  phases  comprised  (1)  a  literature  search  concerning 
particle  formation,  size  change  mechanisms,  and  heterogeneous 
gas/particle  reactions;  (2)  selection  of  size  change  models 
and  three-phase  flow  reacting  gas/particle  coupling  equations; 

C 3 )  coding  of  the  numerical  solution;  (4)  verification  and  docu¬ 
mentation  of  the  code;  and  (5)  preparation  of  a  final  report. 

One-dimensional  nozzle  flow  with  area  change  was  assumed 
in  order  that  the  major  effort  could  be  directed  to  accurately 
modeling  the  physics  of  particle  and  droplet  formation  and  size 
change  mechanisms.  The  overall  accuracy  goal  was  to  obtain 
specific  impulse  predictions  of  condensed  phase  particle  effects 
to  within  ±  1/3%  (1  sec) . 

The  particle  models  were  developed  to  account  for  particle 
agglomeration,  droplet  breakup,  vaporization  and  sublimation, 
melting  and  solidification,  and  condensation,  and  reactions 
that  occur  both  in  the  gas  phase  and  heterogeneously  on  particle 
surfaces . 

Ten  particle  sizes  of  a  single  class  or  different  classes 
are  accommodated  by  the  program.;  Each  class  is  capable  of  inde¬ 
pendent  composition,  phase,  and  size  specification. 

The  reacting  gas/particle  conservation  equations  solved 
account  for  velocity  and  thermal  lags  between  the  droplets/ 
particles  and  the  gas  phase  including  radiative  and  convective 
heat  transfer  to  and  from  the  particles/droplets. 

The  computer  program  calculates  the  entire  nozzle  flow 
field  from  the  nozzle  entrance  to  the  exit  plane.  All  initial 
conditions  are  specified  at  the  nozzle  entrance  including  automatic 
techniques  to  select  the  initial  particle  or  droplet  species 
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and  size  distributions.  The  computation  is  self-contained 
requiring  specification  of  the  propellant  system,  motor  operating 
conditions,  chemical  species  data,  and  nozzle  geometry.  The 
computer  program  is  designed  to  run  with  a  minimum  of  input 
data. 

The  program  has  provision  for  solution  of  frozen  equilib¬ 
rium,  and  nonequilibrium  conditions  as  a  means  of  comparing 
theoretical  specific  impulse  with  predicted  specific  impulse 
losses  due  to  three-phase  flow  effects. 

The  program  will  accommodate  propellant  mixtures  that 
contain  the  following  elements:  carbon,  hydrogen,  oxygen,  nitrogen, 
fluorine,  chlorine,  aluminum,  beryllium,  boron,  magnesium,  zirconium, 
and  iron.  Propellant  formulations  of  primary  interest  to  be 
analyzed  in  detail  have  been  those  that  contain  ammonium  perchlorate 
(NH4C104)  and  aluminum  with  an  N-F  type  binder.  Particular 
attention  has  been  given  to  those  reactions  between  species 
containing  aluminum,  fluorine,  and  oxygen  which  are  predominant 
in  particle  size  change  mechanisms  in  the  formulations  including 
reactions  that  form  A1F,  HF,  OH,  and  A1203. 

The  overall  OD3P  program  consists  of  two  major  programs. 

The  primary  program  is  the  one- dimensional  three-phase  reacting 
gas-particle  kinetic  program.  For  reference  to  theoretical  per¬ 
formance  and  as  a  means  of  obtaining  starting  conditions  for  the 
three-phase  program,  a  one-dimensional  equilibrium  program  has 
been  merged  into  the  three-phase  program.  This  equilibrium 
program  is  designated  herein  as  ODE  (One-Dimensional  Equilibrium)  . 

The  ODE  part  of  0D3P  is  essentially  identical  to  current  NASA  and 
Air  Force  versions  of  this.  ODE  has  been  revised  slightly  to 
write  an  output  file  containing  data  needed  to  start  the  three- 
phase  program. 
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COMPUTER  PROGRAM  STRUCTURE 


To  integrate  the  one-dimensional  gas-particle  equations 
it  is  necessary  to  specify  all  variables  at  the  starting  position, 
furnish  thermodynamic  data  for  all  species  to  be  considered, 
define  the  desired  independent  variable  {temperature,  pressure 
or  area)  as  a  function  of  downstream  position,  specify  integra¬ 
tion  step  and  output  print  controls,  provide  a  set  of  elementary 
kinetic  reactions  for  the  gas  phase  species,  and  specify  the 
number  of  particle  groups,  their  size,  concentration,  composition, 
and  the  size  change  mechanisms  to  be  considered.  The  program 
includes  a  plotting  routine  which  requires  specification  of 
the  variables  to  be  plotted. 

The  major  program  functions  are  arranged  in  seven  modules; 
Figure  1  shows  the  functional  arrangement  of  the  modules  in 
relation  to  the  input  deck  and  several  data  storage  files.  Figure  2 
lists  all  of  the  102  subroutines  in  each  module. 

Table  1  lists  the  input  deck  sections  required  to  run  the 
program.  Table  2  shows  the  format  of  the  input  deck.  Table  3 
lists  the  program  output  parameters  printed  at  print  stations 
specified  by  position  or  area  ratio. 
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INPUT  DECK 


PROGRAM  MODULES 


STORAGE  FILES 


Figure  1.  OD3P  program  structure  showing  relationship  between  input  deck, 
program  modules  and  data  storage  files. 
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Figure  2.  OD3P  overlay  structure 


TABLE  1 

INPUT  DECK  STRUCTURE  AND  FILE  USEAGE 


Input  Deck  Structure  Summary 

Input  to  0D3P  is  divided  into  the  following  sections: 

1.  TITLE  CARD  -  the  contents  of  this  card  will  appear  at  the  top 

of  each  print  station  output. 

2.  $THERM0  -  namelist  input  controlling  input  of  new  thermodynamic 

data  or  the  use  of  a  previously  generated  thermodynamic  data 
file. 

3.  THERMODYNAMIC  DATA  -  includes  gaseous  species  data  and  particle 

species  data  in  either  JANNAF  table  form  or  curve  fits. 

4.  $PROB  -  namelist  for  specifying  operation  of  the  ODE  link 

for  starting  conditions  and  the  ODTP  link  for  the  three-¬ 
phase  flow  calculation. 

5.  REACTANTS  -  ODE  reactant  cards  per  NASA  SP-273 

OMIT  -  omit  species  for  ODE 
INSERT  -  insert  species  for  ODE 
NAMELISTS 

$ODE  -  namelist  input  for  ODE  per  NASA  SP-273  (Ref.  1) 
as  modified  by  SPP  (Ref.  2) 

6.  $CASEIN  -  namelist  data  which  defines  ODTP  initial  conditions, 

the  type  of  calculation  requested,  and  other  specific  case 
input . 

7.  SPECIES  CARDS  -  data  set  containing  the  gaseous  species  to  be 

considered  (in  ODTP)  and  their  initial  concentrations  (if 
not  input  from  ODE) . 

8.  REACTION  CARDS  -  data  set  containing  the  reaction  set  for  the 

gaseous  species,  and  any  third  body  reaction  rate  ratios  to 
be  specified. 

9.  PARTICLES  CARDS  -  when  particles  are  to  be  considered  this 

data  set  specifies  their  species  and  initial  Vpij‘,  Pp^,  Tp*, 
*pk»  Mpk  (note:  p«k  -  lbm  of  particle  per  zt3  of  * 
control  volume,  mPfc  *  bulk  density  ■  lbm  of  particle  per 
ft3  of  particle) .  Note  optional  input  in  $PART  for  ODE 
start. 

10.  $PART  -  when  particles  are  to  be  considered  this  namelist  input 

defines  the  particle  size  change  parameters  specified  for 
the  particles. 

11.  $ADD  -  when  mass,  energy,  or  momentum  addition  for  gas  or 

particles  is  specified  this  namelist  input  defines  the 
schedules  for  the  addition  rates. 

12.  $DATA  -  data  set  for  output  plotting. 

Table  2  shows  the  input  deck  structure. 
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TABLE  2 

0D3P  INPUT  DATA  SET  DESCRIPTION 


Sect  lor. 
Explaining 
Data  Set 


Data  Set 

Cai 
i  : 

i 

Title  card  ! 

rd  Columr 
l)U< 

i 

i  Dumber 
>719 

First  card  » 

so  columns  free  field  title  1 

| 

Last  card 

Won*,  only  i  title  card  allowad  | 

Th 

ermo  Namelist 

| 

1 

First  card 

*  T  H  E 

"MO  | 

Last  card  ”“| 

Send 

1 

CurveTit  Data 

First  cir2 - f 

h  e  r  m 

0  1 

Last  card  E 

K  D 

1 

First  card  l 

OH  T 

**  e  ■  « optional 

Last  card  E 

N  D  L 

ok  T  e  p  h  s  optional  | 

Tabular  Data  1 

First  card 

Last  car4 

Last  particle  H  nr  F  r-mrA  1 

Problem  Card 

First  card 

S  P  R  0 

■ 

Last  card 

1_E.HA 

POE  Reactants  j 

First  card  a 

E  A  C  T 

ANTS 

Last  card  “ 

Blank  card  Required 

ODE  Omit  &  Insert  i 

First  card  n 

MIT 

- 1 

Last  card 

None  Required  j 

First  card  a 

N  S  E  R 

T  1 

Last  card  “ 

None  Reouired  1 

ODE  Namelist 

First  card  n 

A  M  E  L 

I  s  T  S 

Second  card  =”  ~l 

SOD* 

Last  car<k 

SEND 

ODTP  Case  Data 

1  First  card  | 

SC  A  S 

E  1  N  | 

1  Last  card 

SEND 

Species  cards 

First  card  S 

SPEC 

res  1 

- 

Last  card 

none.  Reactions  Must  Follow! 

Reactions 

»  card  p  I 

E  A  r  T 

TONS  1 

I  U«t  card  X. 

AST 

IE. AX  1 

Third  Body  Ratios 

First  card  t 

HIND 

■  n  D  Y  1 

Last  card  i 

AST. 

:  A  A  D  1 

Particle  cards 

1  First  card  p 

A  R  T  1 

:if  s  1 

1  Last  card 

None  Reouired  1 

fra 

rticle  Namelist 

1  First  card 

SPAR 

r 

[  Last  card 

SEND 

Addition  Functions 

1  First  card 

tJ  D  D 

| 

I  Last  card 

SEND 

1 

Plotting  Data 

1  First  card 

S  D  A  T 

i  1 

|  Last  card 

*  EM  P 

_ 1 

2.1 


2.2 


2.2 


2.2 


2.4 


2.5.1 


2.5.2 


2.5.2 


1.5 


2.7 


2.5 


2.5 


2.5 


2.10 

2.11 

2.12 


Notes 


1 


1 


2 


2 


1 


2 


2. 

2. 

4. 

5. 

4. 


Mandatory 

Optional,  not  required  if  data  exists,  select  either  curve  or  table  data. 
Required  if  ODE  option  selected. 

Required  it  ODTP  option  selected. 

Required  if  data  plottinq  selected. 

Optional. 
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TABLE  3 


PRINT  STATION  OUTPUT  DATA 


Gas  Properties 

Velocity 
Temperature 
Heat  Capacity 
Gas  Viscosity 
Density 
Molecular  Wt 
Y 

Prandtl  Number 
Comp ling  Term  A 
Comp ling  Term  B 
Compling  Term  C 

Compling  Term  E 
Mach  Number 
Pressure 
Area 

Axial  Position 
Time 

Iterations 

Delivered  Thrust 
Coefficient 

Delivered  Spec. 
Impulse 

Vacuum  Thrust 
Coefficient 

Vacuum  Spec. 
Impulse 

C*t  Characteristic 
Velocity 


ft/sec 

•R 


BTU/lbm-*R 

lbm/ft-sec 

lbm/ft3 


lbm/lbmole  (of  gas  only) 
unitless 


unitless 
1/Unit  Length''* 

1/Unit  Imngth|see  Analysis 

lbm/ (f t-sec  - /for  definition  of  the 

Unit  Length)  [  coupling  terms. 

1/Unit  Length) 
unitless 


PSIA 

ft2 

Unit  Length 
seconds 
unitless 
uni/tless 

sec 

unitless 
sec 

ft /sec 


Downstream  of 
Throat  only 
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TABLE  3  (contd) 


Species  Concentra¬ 
tion 

Gass  Mass  Flux 

Particle  Mass  Flux 

Total  Mass  Flux 

Cumul.  Mass  Flux 
Error 

Cumul.  Mass  Flux 
Error 

Gas  Energy  Flux 

Particle  Energy 
Flux 

Total  Energy  Flux 

Cumul.  Energy  Flux 
Error 

Cumul.  Energy  Flux 
Error 

Gas  Momentum  Flux 

Particle  Moraen.  Flux 

Total  Momen.  Flux 

Cumul.  Momen.  Flux 
Error 

Cumul.  Momen.  Flux 
Error 

Current  Step  Size 

P  used/p  -  RO.R.T 

Gas  H  (H  +  V2/2) 

Molecular  Weight 

Molecular  Weight 

Integrated  Gas 
Mass  Added- 

Integrated  Part. 

Mass  Added 

Integ .  Energy 
Added 


Mass  Fraction 

lb/ sec 
lb/ sec 
lb/sec 
lb/sec 

percent 

Btu/sec 

Btu/sec 

Btu/sec 

Btu/sec 

percent 

ft  -  lbra/sec2 
ft  -  lbm/sec 

y 

ft  -  lbm/sec 
ft  -  lbm/sec2 

percent 

unit  length 

uni,tless 

Btu/lbro 

lb  mixture/mole  gas 
lb  mixture/mole  mixture 
lbm/sec 

lbm/sec 

Btu/sec 


TABLE  3  (contd) 


Wall  Radiation 
Flux 

Wall  Temperature 

Integ.  Momentum 
Added 

1  -  Ic.,  Species 
Summation 

Percent  dT/T 
per  step 

Energy  Error  per 
step  as  dT/T 

Maximum  Species 
dci/ci 

A. dx 

B. dx 

Gas  Equivalence 
Ratio 

Integ.  Particle 
Mass  Transfer 

Mean  Diameter , 
all  particles 

Geometric  Standard 
Deviation,  all 
groups 

Mean  Diameter, 
each  set 

Geometric  Standard 
Deviation,  each 
set 


Btu/s/ft2  wall 
•R 

ft/sec 

mass  fraction 
#R/®  R  per  step 
percent  per  step 
percent  per  step 

unitless 

unitlesB 

unitless 

lb/s/ft3  gas 

micrometers 

unitless 

micrometers 

unitless 
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TABLE  3  (contd ) 


Particle  Properties  (Printed  for  each  Group) 
Particle  Phase 


Radius 

ft.,  microns 

Phase  Volume 

Fraction 

unitless 

Velocity 

ft/sec 

Density 

lbm/ft^  volume  of 

Temperature 

•R 

Bulk  Density 

lbm/ft3  particle 

Wpe,  Wp1,  wpc,  Wpk 

lbm/(ft3  sec) 

Number  Density 

Partieles/ft3 

Compling  Term  D 

1/unit  length 

Compling  Term  F 

1/unit  length 

Reynolds  Number 

unitless 

Nussult  Number 

unitless 

Knud sen  Number 

unitless 

Heber  Number 

unitless 

Drag  Coefficient 

unitless 

Velocity  Relaxation 

(sec-Unit  L)/ft 

Temperature 

Relaxation 

(sec-Unit  L)/ft 
• 

V  •  V  A 

lbm/sec 

Particle  Mass  Flux 
Error  Term 

lbm/sec 

Integrated  Particle 
Vaporization 

lbm 
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PARTICLE  BREAKUP  AMD  COLLISIONS 

Particle  breakup  calculations  using  the  Extended  Delta 
SPP  test  case  were  performed  starting  with  30  pm  and  300  pm 
median  particle  sizes.  Nine  groups  were  used  with  an  initial 
geometric  standard  deviation  of  1.9.  Table  4  shows  that  for 
both  calculations  the  particles  broke  up  to  nearly  the  same 
mean  diameter.  Of  interest  is  that  the  larger  initial  sizing 
(300  pm)  produced  a  slightly  smaller  final  diameter.  The  breakup 
ratio  used  was  2  to  1,  that  is,  at  any  point  where  a  particle 
Weber  number  exceeds  the  critical  value  of  4/CD,  the  particle 
diameter  is  reduced  by  a  factor  of  2.  Figure  3  shows  the 
size  distributions  for  the  initial  30  ym  median  size  case  and 
at  the  throat  for  that  case.  The  throat  distribution  remains 
approximately  log-normal  as  expected  but  the  geometric  standard 
deviation  is  much  lower.  There  is  also  a  shift  in  the  size 
arrangement  of  the  groups.  Groups  1,  5,  and  9  remain  respectively 
the  smallest,  median  and  largest  size  but  other  groups  shift 
around.  The  geometric  standard  deviation  of  1.9  for  the  initial 
particle  sizing  is  actually  based  on  measurements  at  the  nozzle 
exit.  It  would  appear  that  a  larger  value  should  be  used  at 
the  nozle  inlet  so  that  the  exit  distributions  would  be  closer 
to  1.9. 

Figure  4  taken  from  NASA  SP8039,  shows  particle  size 
correlations  for  solid  rocket  motors.  The  Extended  Delta  test 
case  mean  particle  sizes  agree  well  with  the  NASA  SP8039  curve 
but  is  higher  than  the  SPP  correlation  (which  is  currently  being 
revised) . 

Particle  size  as  a  function  of  axial  position,  plotted 
by  OD3P,  is  shown  in  Figure  5  for  an  initial  D^j  of  30  pm 
with  nine  particle  groups  and  a  geometric  standard  deviation 
forsize  distribution  of  1.9.  The  position  X  ■  0  is  the  throat. 
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TABLE  4.  EXTENDED  DELTA  TEST  CASE  PARTICLE  BREAXUP 

AND  COLLISIONS 


Case  1 

Case  2 

Initial  median  diam,  ym 

30 

300 

Initial  mean  diam,  D^,  ya 

36 

360 

Initial  geom.  std  dev. 

1.9 

1.9 

Throat  mean  diam,  D^,  ym 

8.45 

7.51 

Throat  geom.  std  dev. 

1.26 

1.232 

Exit  mean  diam,  D^3,  ym 

8.34 

7.14 

Exit  geom.  std  dev. 

1.32 

1.288 
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A  Ext.  Delta,  30  pm  by  OD3P 
vs.  D.  □  Ext.  Delta,  300  ym  by  OD3P 


Figure  4.  Particle  size  comparison  for  breakup  calculation 


All  nine  groups  break  up  before  or  at  the  throat;  some  groups 
break  up  more  than  once.  Particle  size  of  the  larger  particles 
beyond  the  throat  grows  slightly  due  to  collisions.  Collision 
rate  is  low  because  the  particles  are  all  nearly  the  same  size. 

SPECIES  AND  KINETIC  REACTIONS 

KVB  developed  a  species  and  kinetic  reaction  set  for  a 
typical  fluorinated  propellant  containing  23  gaseous  species  and 
76  reactions.  This  set  was  presented  at  the  12th  JANNAF  PSS  meeting 
in  Salt  Lake.1  Since  that  meeting  the  set  of  76  reactions  has  been 
updated,  as  shown  in  T:-.ble  5,  to  include  the  most  recent  kinetic 
rate  data. 

In  using  this  set  it  is  important  to  understand  the 
assumptions  on  which  it  is  based.  The  reaction  set  was  developed 
for  a  specific  set  of  species  as  shown  in  Table  5.  Species  are 
shown  in  four  groups.  The  first  group  are  those  that  are  important 
to  the  energy  release  during  nozzle  expansion  of  a  non-fluorinated 
propellant.  These  species  are  listed  in  order  of  decreasing 
magnitude  of  energy  release.  The  second  group  shows  additional 
major  and  minor  species  for  fluorinated  propellants.  The  third 
group  includes  species  that  are  moderately  high  in  concentration 
(>0.0001  mole  fraction)  but  do  not  contribute  significantly  to 
energy  release.  The  fourth  group  contains  three  species  which 
are  not  important  to  energy  release  but  are  important  to  consider 
for  performance  based  on  restricted  equilibrium.  This  is  due  to 
the  fact  that  for  full  equilibrium  AlCl^  does  not  form  during 
expansion;  rather  it  disappears.  However,  when  AljO^  formation 
is  suppressed  as  in  restricted  equilibrium,  fairly  large  amounts 
of  AlCl^  are  formed. 


Cherry,  S.,  "Kinetic  Rate  Data  Screening  and  Update  Procedure," 
12th  Meeting  Minutes  of  the  JANNAF  Performance  Standardization 
Subcommittee,  Salt  Lake  City,  Utah,  January  25-26,  1979. 
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In  developing  the  set  of  76  reactions,  the  following 

23  species  were  used:  15  species  from  the  first  group  (excluding 

Al2<>3  and  N2) »  the  five  major  fluorine  species  from  the  second 

group,  and  N  and  A10H  from  the  third  group.  The  latter  two  were 

included  to  provide  a  kinetic  path  for  NO.  Selection  of  these 

species  was  based  on  ODE  calculations  only  with  fluorinated 

propellants.  Therefore,  the  use  of  a  subset  of  the  76  reaction 

set  for  non-fluorinated  propellants  requires  further  consideration. 

At  the  12th  PSS  meeting,  Steve  Cherry  of  KVB  presented  a  preliminary 

screened  set  of  16  reactions  involving  18  gaseous  species,  excluding 

fluorinated  species.  This  set  was  based  on  previous  screening 

of  an  original  Al  set  included  in  the  SPP  initial  version.  A 

further  subset  of  fifteen  reactions  (excluding  CO  +  0  *  02)  with 

18  species  was  screened  using  updated  rates  from  the  final  KVB  set 

2 

by  Software  Engineering  Associates  in  the  SPP  program.  These 
results  indicated  a  change  in  restricted  equilibrium  Ifip  of 
only  0.08%  compared  with  the  full  set  of  reactions  shown  in 
Table  5  (excluding  fluorinated  species  and  including  four  reactions 
involving  Cl2,  A1CI»2,  and  AlCl^).  This  SEA  set  of  15  reactions  is 
therefore  considered  appropriate  for  use  with  non-fluorinated 
propellants  when  performing  restricted  equilibrium  calculations. 

The  fact  that  this  set  allows  AlCl-j  to  increase  during  expansion 
requires  further  assessment  when  the  restriction  on  mass  transfer 
is  removed  as  in  0D3P.  KVB  plans  further  work,  starting  from  the 
SEA/KVB  set  of  15  reactions  and  excluding  the  reactions  involving 
Cl2,  A1C12,  and  A1C13. 


2 - 

Nickerson,  G.  R. ,  Coats,  D.  E.  and  Hersen,  R.  W. ,  "Solid  Rocket 
Motor  Performance  Predictions  Using  the  Improved  SPP  Computer 
Model,"  16th  JANNAF  Combustion  Meeting,  Monterey,  California, 
September  10-14,  1979. 
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TABLE  5.  REACTION  SET  FOR  ALUMINIZED 
PROPELLANTS  WITH  N-F  BINDERS 


REACTIONS 

N,N«N2, 

Aat.ElB, 

BaO.O, 

CHERRVU967) 

001 

H»H»H2* 

A»l,09tl», 

Hal,, 

*■0,0, 

JEN3En/JUNES(197B) 

00? 

H,CL»mCl* 

Aal  ,<I5E  22* 

►•■2., 

B»0.O, 

JEn3En/jOnE*(197*> 

001 

ALOF»*LF,D, 

A»9,5E10, 

*■•0,3, 

*■110,0, 

B  AHN (1970) 

00a 

0*H»0h# 

ia3,»2ElM, 

*■!., 

0*0,0, 

JEnSEn/ JO*E 5(1 978) 

005 

0H*m«h20» 

•■3.22E22, 

*■*., 

*■0,0, 

JEnSEn/ JOnE 3(1978) 

Otlb 

*LCuF»*U^»Ct» 

AaA.SElO, 

Na.0,5, 

B»B 1,2, 

8AHN( 1970) 

007 

Aa| ,6E12, 

*■•0.3, 

B«113. 3, 

*AHN(1970) 

00* 

N,0»NU» 

Aaa.aEtb, 

*■0.9, 

BaO.O, 

BaULCmI 19731 

009 

AL0»*L*O» 

•■1 ,7E 12, 

*■•0,5, 

*■113.3, 

*ANN(|970) 

010 

AL0*CH«itP2M. 

•■3,0Ei6, 

*«0.3, 

*■0,0, 

estimate 

oil 

ALOCL«*LO*Cl> 

•■O,bEl0, 

*■•0,5, 

*■120,3, 

BamN( 1970) 

•  12 

AL0CL>AlCL«O. 

AaO.bElO, 

*■•0,3, 

*■120.3, 

*AHN(1970) 

013 

CO«Q*C02» 

Ab2.5<*E  15, 

*■0,0, 

*•0,37, 

JEnSEn/ JUNE 5 ( 197* ) 

014 

alo*n«alon. 

•■3,0E lb. 

*■0,5, 

*■0,0, 

estimate 

•  15 

AL*OH»*LON, 

•■3.0E16, 

*■0,5, 

*■0,0, 

ebtinate 

016 

End  T8R  REAX 
CU*H2»NCU»H, 

AaB.OSEl?, 

*■0.0, 

B»0,233, 

JEN*EN/ JONES (1978) 

017 

H260»0M*M, 

Aal.BOElO, 

*■>1,0, 

*■8.903, 

JEn5En/J0*E5(197B) 

01* 

M2*0N»N20*Nf 

A»1.1«E9, 

*■-1.3, 

*■3.627, 

JEnSEn/ JONES (1978) 

019 

C0*0"»C02*H, 

Aal.b9E7, 

*■•1 .3, 

*■-0.656, 

JEn9En/j0nES(197B) 

020 

CL*0h«hCL*0, 

•■2,«lE12, 

*aO,Q, 

B»0,96B, 

JENSEN /JONES (1978) 

021 

HCL60H»M20*CU, 

•■1.3QE13, 

*•0,0, 

*■2,0*7, 

JEN9En/J0nES(197B) 

022 

0H40H»0tH20, 

•■6,30112, 

*■0,0, 

*■1,093, 

JENSEN/ JONES (1978) 

023 

C02*N2»C0*Hj0# 

IM.SU, 

*■-0,3, 

B»iS„ 

TUNDER(l9b7) 

02o 

N»0*»N0»M, 

•■S,3Elt, 

*■-0,3. 

*■5,628, 

Cherry (i9o7) 

025 

K2*OaNO,N, 

*■7, 60113, 

NaO.O, 

*■75,516, 

JENSEn/JUNES( 1970) 

02* 

C02*N»C0,N0, 

•■1.05EI1, 

*■-0,9, 

*■59,62, 

CHERRY! 1967) 

027 

A' -0“-*'  n*M, 

*■9. Ell, 

*■-0,5, 

*■5,619, 

ESTIMATE 

028 

AlO»CO»AL*Co*» 

•■l.Ell, 

*■-0,3, 

•■6.066, 

estimate 

029 

al*mcu*alCl*h, 

•■5. Ell, 

*■•0,5, 

*■5.673, 

estimate 

•30 

Ac0*CLRAlCl40, 

•■9,111, 

*••0,5, 

*■6,066, 

estimate 

•31 

ALCL»C02»*L0C|.*C0f 

•■1,111, 

*»-0;5» 

*■0,995, 

estimate 

•32 

ALO»HCLrAU)cL*H, 

•■l.Ell, 

*■•0,5, 

•■5,673, 

estimate 

•33 

ACCL»OH«ALOct*M, 

•■l.lll. 

*■-0,3, 

•■5,619, 

Estimate 

•3a 

MFAAcRAtfAH. 

•■9, Ell, 

N»-0,3, 

*■7,006, 

estimate 

•33 

AuE2»M*Atr»HE, 

•■9,111, 

*■•0.3, 

■■7,175, 

estimate 

•36 

AUCLA*MbA17,MCC, 

•■9, El 1 , 

*■•0,3, 

*■0,707, 

estimate 

•37 

ALCL^*H»H^tALCt» 

••3, Ell, 

*»-0,5» 

*■0,961, 

estimate 

•38 

(continued) 
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TABLE  5  (Continued) 


ALCtE*0«*L0f*CL» 

AaS.Cll  , 

Na-0,5, 

BbA.TOT, 

CSUhATC 

•so 

Om**LOcL,CL**LOS**» 

AbI.EU, 

NB-0.5, 

SbT.ITO, 

estimate 

060 

N«ALObnO«AL, 

*■5,111. 

Nb-0.5, 

1'i.m, 

cstxhatc 

061 

ALO»AiCl"Ai*AIOCi# 

A*1 ,E 11 , 

Na-0,5, 

Dak, SIS, 

estimate 

••2 

C02«*L^*CO»ALOIr# 

Aal.Ell, 

Nb-0,5, 

SBk.OOS, 

estimate 

•B| 

ALE«HCL«HF«*LCLt 

AbI.CIO, 

Nb-O.S, 

BaTS.TSS, 

estimate 

•BA 

NT  4 ALOa ALUP+H* 

AiliUlf 

Na-0.5* 

Ba7. Aik, 

EBTIMATC 

•«5 

hE«alubalt«on» 

ABl.ClO, 

Nt-0,5, 

0-71,025, 

estimate 

06k 

Al02*hClbHT4AlCl0, 

ABl.ClO, 

Nb-O.S, 

BBkS.AIt, 

estimate 

•67 

HTAAi.uCL*ALor«MCu 

AbI.CIO, 

Na-O.S, 

Bb7S,I21, 

estimate 

•as 

AlO*nCl*OH*ALCL# 

AB1.E10, 

NB-O.S, 

Bakt.700, 

estimate 

•so 

ALAALF2«ALr4ALE, 

AaS.Cll, 

Na-O.S, 

BbT.175. 

Estimate 

•so 

AlT»OX»AIOT4H, 

AbI.EU, 

Na-O.S, 

SaS.klO, 

estimate 

OSl 

alClfbalbalfbalCl# 

abs.cii. 

Na-O.S, 

BBk.Okl, 

estimate 

•52 

ALT«ALOaALOT«ALf 

AbI.EU* 

Na-O.S, 

•■7**60, 

estimate 

•SI 

ALT2*ALO»ALT*AlOT, 

AbI.EU, 

Na-O.S, 

BbT.175, 

ESTIMATE 

•SB 

AL^«AlClEbAlCLbAI.F2» 

ABliUll 

Na-O.S, 

BBk.Okl, 

estimate 

OSS 

ALCLFtALOBAL^aALOCL# 

abi.eu. 

NB-O.f, 

Baa.707, 

estimate 

•5k 

ALT  a ALOCLb AlCLb ALOT  , 

AbI.EU, 

Na-O.S, 

•■7.175, 

estimate 

057 

ALOCL*RBALCtARO« 

abS.CU, 

Na-O.S, 

BB7.17S, 

estimate 

OSS 

ALCLF  t ALOa ALCLBALOT  , 

ABl.lUf 

Na-O.S, 

Bak.Okt, 

ESTIMATE 

•so 

AL0*H20>AL02MtM( 

ABi.cn, 

NS-O.S, 

BBk.Skk, 

Estimate 

OkO 

ALObnbnOaAL, 

abs.cu, 

nb-0,5. 

B-k.Akk, 

estimate 

•61 

alF2*alocl*alclf*alof, 

Aal.ElO, 

Nb-O.S, 

0*73.050, 

Estimate 

0k2 

AL0F4N«ALT«n0» 

A«s,cn, 

nb-0,5. 

8*7,500, 

estimate 

OkS 

AL04NCL«AL0H«CL» 

Aai.cn, 

Nb-O.S, 

B*S,k7S, 

ESTIMATE 

•ka 

ALCL«0M"AL0H4CL* 

ABI.EU, 

Na-O.S, 

•■k.Skl, 

estimate 

•65 

AL0H40M»N*Ai.02M, 

abi.eh. 

Na-O.S, 

6*5,027, 

estimate 

•66 

ALAM2oaALOH«Mt 

abs.ch. 

Nb-0,S, 

BBk.Skk, 

ESTIMATE 

067 

AL0«H2«AL0M4H, 

ABi.cn, 

Nb-O.S, 

0-5. 731, 

estimate 

OkS 

ALO«OMbALOM40> 

abi.cii* 

Nb-O.S, 

B«S. 627, 

ESTIMATE 

060 

AL0H4NP«aL0P*H2» 

AbI.CIO, 

Na-O.S, 

•■70,33k, 

estimate 

•70 

AL0M4HP»Atr,H20, 

ABi.cie, 

NB-O.S, 

•*7B,k*5, 

ESTIMATE 

071 

AL0"*«CloAL0CL«H2( 

ABl.ClO, 

Nb-O.S, 

•■ki.no. 

Estimate 

•72 

AL0H4HCL»ALCL*m20, 

ABl.no, 

Nb-O.S, 

•BkS.OlO, 

estimate 

07S 

ALP*AL02N»AL0H4At0F# 

ABI.CIO, 

NB-O.S, 

•■•3.030, 

ESTIMATE 

07a 

ALON4H20«AloXM4H2( 

ABl.ClO, 

NB-O.S, 

•■•S.657, 

estimate 

07S 

AL02N4A|.CLBAL0N«AL0CL» 

ABI.CIO, 

nb*O.S, 

•■72.656, 

estimate 

07k 

i«IT  RCtl 

fNlRO  SOOT  R[U  RATE  RATIOS 
ALL  EQUAL  |,0 
LAST  CARD 
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TABLE  6.  SPECIES  CONSIDERED  IMPORTANT  FOR 
SOCKET  MOTOR  PERFORMANCE  ANALYSIS 


X.  MAJOR  SPECIES  BASED  ON  ENERGY  RELEASE  (NON-FLUORINATED) 


A12°3 

HCI 

0 

A10C1 

CO 

OH 

A1C1 

O 

O 

M 

N2 

H2° 

NO 

H 

Cl 

AlO 

H2 

A1 

aio2h 

ADDITIONAL 

SPECIES  IMPORTANT  FOR 

FLUORINATED  PROPELLANTS 

4 

Mai  or 

Minor 

HF 

aicif2 

aif2 

F 

AlOF 

aici2f 

A1F 

aif3 

A1C1F 

III.  SPECIES  HIGH  IN  CONCENTRATION  BUT  NOT  IMPORTANT  TO  ENERGY 
RELEASE 


A1H 

°2 

A1C12 

N 

ai2o 

AlOH 

IV.  SPECIES  IMPORTANT  IN  RESTRICTED  EQUILIBRIUM  BUT  NOT  FOR 
FULL  EQUILIBRIUM 


aici3  aici2  ci2 
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CONSTANT  FRACTIONAL  LAG  PROCEDURE 
FOR  DRIVING  VARIABLE 

One-dimensional  kinetic  calculations  required  specification 
of  a  defined  independent  variable,  usually  pressure  in  the  sub¬ 
sonic  and  sonic  regions  and  area  in  the  supersonic  region.  The 
pressure  function  is  generated  using  an  average  equilibrium  pres¬ 
sure  expansion  coefficient  from  the  chamber  to  the  throat.  This 
is  an  appropriate  approximation  for  calculations  with  gas-particle 
velocity  and  temperature  equilibrium.  However,  for  the  non¬ 
equilibrium  procedure  of  003P,  it  was  necessary  to  revise  the 
calculation  to  include  an  approximation  more  appropriate  for 
flow  with  particle  velocity  and  temperature  lags  relative  to  the 

gas.  The  approximation  used  is  the  constant  fractional  lag  pro- 

3 

cedure  of  Kliegel  .  The  procedure  was  developed,  as  shown  in 
an  appendix  to  this  paper,  for  both  pressure  and  temperature 
defined  calculations.  While  the  pressure  defined  mode  is  con¬ 
ventional,  the  temperature  defined  mode  is  novel  and  was  developed 
to  improve  the  numerical  solution  of  gas-particle  flows  with 
kinetics.  Solution  of  the  gas-particle  equations  in  OD3P  involves 
up  to  104  simultaneous  equations  for  the  gas  species  concentrations 
(40);  gas  velocity,  density,  area  and  pressure  or  temperature  (4); 
and  six  equations  for  each  of  10  particle  groups  (60) .  For  near¬ 
equilibrium  flow,  numerical  instabilityis  a  major  difficulty 
in  obtaining  solutions  because  the  species  kinetic  rates  are 
such  strong  functions  of  temperature.  The  use  of  a  temperature- 
defined  mode  for  specifying  the  driving  variable  was  found  to 
minimize  any  tendency  for  numerical  instability  due  to  the  stiff 


Kliegel,  J.  R. ,  "Gas  Particle  Nozzle  Flows,"  Ninth  (International) 
Symposium  on  Combustion,  The  Combustion  Institute,  1962. 
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equation  problem.  This  procedure  is  considered  useful  only  for 
cases  in  which  the  gas  mixture  is  in  fact  very  near  equilibrium 
and  the  effect  on  prediction  of  kinetic  loss  has  not  yet  been 
evaluated  in  full  gas-particle  calculations.  For  cases  where 
significant  kinetic  non-equilibrium  exists,  it  is  considered 
more  appropriate  to  retain  the  pressure-defined  mode. 

ROCKET  MOTOR  VERIFICATION  CASES 

In  order  to  assess  the  accuracy  of  0D3P  and  demonstrate 
its  operation,  a  set  of  six  motor  test  cases  is  being  used. 

These  cases  are  identical  to  those  currently  being  used  for  the 
Improved  SPP.  Final  execution  of  these  cases  is  not  yet  completed 
pending  accuracy  verification  of  all  the  particle  size  change 
models  in  0D3P,  together  with  kinetics  and  numerical  integration 
methods.  However,  all  cases  have  been  run  successfully  through 
the  ODE  module  and  through  the  initial  setup  routines  of  OD3P, 
including  generation  of  tables  for  either  pressure  or  temperature 
defined  independent  variable.  Table  7  presents  the  lag  parameters 
that  result  from  the  independent  variable  generation.  Particle- 
to-gas  mass  flow  ratios  (frD/Sr  )  are  similar  for  all  the  motors 
except  RSM  which  is  very  low.  The  nozzle  discharge  coefficients 
are  also  similar  for  the  first  five  cases  but  lower  as  expected 
fcr  P.SM.  All  the  results  shown  were  based  on  particle  size  using 
the  Improved  SPP  particle  size  correlation.  The  particle  size 
from  that  correlation  is  checked',  against  a  particle  breakup 
criteria  using  a  Weber  number  of  28.  In  all  cases  the  particle 
size  from  the  correlation  was  below  the  critical  breakup  size. 

The  prediction  of  CD  needs  to  be  further  assessed  based  on  actual 
0D3P  particle  calculations  which  are  yet  to  be  completed. 
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TABLE  7.  CONSTANT  FRACTIONAL  LAG  PARAMETERS 


Motor 

£ 

1 

Y 

K* 

L* 

CD 

Extended  delta 

0.3979 

1.1247 

0.886 

0.871 

1.0150 

Titan  III  C 

0.4116 

1.1288 

0.843 

0.767 

1.0240 

C4  Stage  3(ADP) 

0.4669 

1.1237 

0.839 

0.767 

1.0276 

IUS  Large  Motor 

0.4708 

1.1194 

0.834 

0.767 

1.0274 

AIM 

0.4264 

1.1255 

0.835 

0.753 

1.0257 

Reduced  Smoke 

0.004 

1.1675 

0.821 

0.683 

1.0031 

Maverick 
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APPENDIX  A 


CONSTANT  FRACTIONAL  LAG  METHOD 
FOR  TEMPERATURE  AND  PRESSURE 
DEFINED  SOLUTIONS 

Obtain  the  following  from  equilibrium  solution: 
Tc  chamber  temperature 
T* ,  throat  temperature 
h_  chamber  enthalpy 

C  t 

h* ,  throat  enthalpy 


Calculate  equilibrium  expansion  coefficient. 
Tc 

Ye  =  2  ^  -  1 


Calculate  average  specific  heat 

c  .K^l 


T  -T 
C 


Calculate  average  gas  specific  heat 


C  =  C 

pg  p 


(■  %Y 


pp 


(C-l) 


(C-2) 


CC-3) 
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gjgij 


T 


Calculate  gas  expansion  coefficient 


(C-4) 


Calculate  the  two-phase  flow  expansion  coefficient  by  iterating 
the  following: 

/•  b  \*  *2 


C=l+ 


Y  »  1+ 


x 

>  f A7(w^F 

"  (“p/ws)*  !k*  [(l-K‘)  VK‘]  +  (v1)  (cpf/^pg  )*“•*}  ,c-6' 


(V1)* 


(C-7) 


In  the  above,  assume  the  particle  lags  are  zero  to  begin  the 
iteration;  i.e.,  K*«=  L*  =  1 

Estimate  the  two-phase  flow  throat  conditions 


*  T 

*g  -  2  9 

T  +  1 


\  B  Y  +  1/ 

•A)* 


CC-8) 


CC-9) 


(C-10) 
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For  cases  with  no  particles  the  above  equations  are  sufficient  and 
no  iteration  is  necessary. 

For  flow  with  particles  the  following  equations  provide  an  estimate 
of  the  particle  lags  at  the  throat. 

Calculate  the  mass  average  particle  size  for  particles  considered 
in  the  case 

_  k  W/" pk  (C-ll) 

43  k  Opk/V 


Calculate  a  maximum  particle  size  based  on  particle  break  up 


W 

*  gc  °  ecrit 
Crit  2pg  Vg  <1-K*> 


(C-12) 


Set  rp*  to  the  minimum  of  r43  and  r 

Calculate  CD*,  fp*,  y*  and  Pr*  by  standard  formulas. 
Calculate  particle  density  at  the  throat  from: 


(C-13) 

(C-14) 

(015) 

(016) 
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Calculate  the  estimated  particle  lags. 

*  *  * 

P  fp  r 
D*  s  * — fc- 


*  l  ¥ 

m  r  V 
P  P  9 


(C-17) 


*  Vp*  9  * 

K  a  -Ey  =  J  D 

IT  ” 

2  r* 

1/2 

[-A 

L  9D 

(¥  * 

1/2 

-1 

vg 

r 

X  '  J 

(C-18) 


L*  *  rJ1  *  [l  +  3  Pr*  (cpp/£pg  )  *(1~K*)/K*]"1  ,c-19’ 

c  g 


Iterate  the  entire  solution  procedure  above  until  K*  and  L* 
converge. 

After  convergence,  calculate  predicted  nozzle  mass  flow 


-  [3  +(“p/“ a)*]  *  r‘2  V 

Calculate  the  inlet  contraction  ratio 


(C-20) 


(C-21) 
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Calculate  a  temperature  table  as  a  function  of  nozzle  position  with 

N  values  between  the  inlet  and  the  throat. 

* 


AT 


Table 


(C-22) 


Calculate  the  area  ratio  for  each  table  value  of  T,. 


a 

-  .  T 

Y-l  g 

_i_  !c 

Y-l 

“ «■  - 
a  i 

2  T  -T 
c  g 

Y+l  Tg 

1/2 


(C-23) 


Calculate  a  table  of  nozzle  position  from  a/a*  and  the  nozzle 
geometry.  This  produces  a  table  of  T^  and  a  versus  x.  The 
table  of  Tg  and  a  versus  x  can  be  used  in  two  ways  during  the 

temperature-defined  solution  at  the  two-phase  flow: 

1.  Use  Tg  =  f  (x) ,  where  x  is  the  local  position  during 
solution. 

2.  Use  T  *=  f  (a) ,  where  a  is  the  local  calculated  area . 

The  first  method  is  the  most  stable  because  Tg  will  remain  fixed 
at  the  forward  point.  However  there  will  be  a  difference  between 
calculated  area  and  the  tabular  area.  The  throat  may  occur  at  a 
point  that  does  not  match  the  Tg*  calculated  above. 

The  second  procedure  wil;l  tend  to  be  less  stable  because 
the  forward  point  temperature  will  change  as  a  function  of  the 
calculated  area.  However,  the  calculated  area  may  better  match 
the  geometric  area. 

The  first  method  was  selected  for  the  program. 
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The  above  procedure  is  used  to  generate  a  temperature 
table  for  a  temperature  defined  solution  which  is  very  stable 
for  numerical  integration  when  gas  phase  chemistry  is  near 
equilibrium  and  the  gas  species  differential  equations  are  very 
stiff. 

The  temperature  defined  procedure  can  be  converted  to 
pressure  defined  by  the  use  of 


(Y/Y-l) 


(024) 


and 


dP 

dx 


\  _  Y  T 

Y/(Y-1>  -  1 
c  ' 


dT 

dx 


(025) 
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TABLE  6.  SPECIES  CONSIDERED  IMPORTANT 
ROCKET  MOTOR  PERFORMANCE  ANALYSIS 
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EFFECT  OF  MASS  TRANSFER  ON  PERFORMANCE 
EXTENDED  DELTA  TEST  CASE 


CONSTANT  FRACTIONAL  LAG  PARAMETERS 


Motor 

% 

Y 

K* 

L* 

CD 

Extended  delta 

0.3979 

1.1247 

0.886 

0.871 

1.0150 

Titan  III  C 

0.4116 

1.1288 

0.843 

0.767 

1.0240 

C4  Stage  3(ADP) 

0.4669 

1.1237 

0.839 

0.767 

1.0276 

IUS  Large  Motor 

0.4708 

1.1194 

0.834 

0.767 

1.0274 

AIM 

0.4264 

1.1255 

0.835 

0.753 

1.0257 

Reduced  Smoke 

0.004 

1.1675 

0.821 

0.683 

1.0031 

Maverick 


OD3P  RESULTS  FOR  EXTENDED  DELTA 


Inert  particles 
0*3-3.341 


Mass  transfer  only 
D*3=3.33 


All  size  change 

* 

D^ 3=6. 26pm 


1  rVro 

K* 

L* 

1.55 

.962 

.910 

2.86 

.911 

.805 

5.33 

.832 

.657 

1.58 

.962 

.879 

2.88 

.911 

.773 

5.34 

.833 

.622 

1.53 

.960 

.902 

3.19 

.895 

.771 

7.03 

.796 

.564 

A  RtMfch-Cottrtl  Company 
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Appendix  12:  Particle  Impact  Erosion  Program 
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PERFORMANCE  PREDICTION  METHODOLOGY 
STATUS  SUMMARY 


! 

|  Joe  D.  Hoffman 

Purdue  University 


February  1980 


Performance  Prediction  Methodology  Status  Summary 
Performance  -  Specific  impulse,  Isp 

Theoretical  performance  -  Maximum  Isp  with  no  losses,  Isp^h 
Losses  -  Decrement  in  Isp  due  to  real  effects,  Alsploss 
Delivered  performance  —  Isp  including  losses,  IspD 

The  objective  of  the  performance  prediction  methodology  is  to  a  priori 
predict  IspD  within  ±  0.5%. 

Theoretical  performance,  Ispth>  is  determined  by  a  one-dimensional,  isentropic 
(i.e.,  adiabatic  and  frictionless),  equilibrium  expansion  from  equilibrium 
combustor  conditions.  This  calculation  is  accomplished  with  the  NASA 
Lewis  thermochemistry  program  known  as  ODE. 


LOSSES 


Two-dimensional  loss 

0.1  to  3% 

Finite  rate  kinetics  loss 

0.05  to  10% 

Two-phase  flow  loss 

0  to  5% 

Boundary  layer  loss 

0.5  to  5% 

Erosion  loss 

0  to  2% 

Submergence  loss 

0  to  1% 

Combustor  efficiency 

0  to  5% 

lsp„  ■  lspth  -  l  41spl0ss 

IspD  =  lspth  v nlQSS 

IspD  -  lspth  -  HspBL 
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Loss  Calculations 

Loss  Model  -  Performance  calculation  with  a  model  containing  certain 
mechanisms  including  a  particular  loss. 

Reference  Model  -  Performance  calculation  with  a  model  containing 
the  same  mechanisms  except  the  loss  in  question. 

Problems  -  Selection  of  the  losses  to  be  separated  or  coupled. 

Establishment  of  the  loss  models  -  very  subjective. 

How  to  make  references  model  calculations. 

I 

I 


COMBUSTION  EFFICIENCY 


KINETIC  EFFICIENCY 


Advanced  Transonic  Program 

The  Advanced  TranSonic  (ATS)  program  determines  the  subsonic  and 
transonic  flow  field  in  a  solid  propellant  rocket  motor. 

Unsteady  two-dimensional  flow 
Steady  flow  solution  at  large  time 
Coupled  gas-particle  flow 
Perfect  gas 

Constant  size  particles 

The  result  of  the  ATS  program  is  a  supersonic  initial-value  downstream 
of  the  nozzle  geometric  throat  from  which  the  TD2P  analysis  can  be 
initiated.  A  portion  of  the  nozzle  submergence  loss  is  included  in 
the  results  predicted  by  the  ATS  program. 


New  Considerations 


One-Oimensional  Three-Phase  Flow  Program 

0D3P  -  particle  radiation 

solidification  and  crystallization 
collisions  and  agglomeration 
fragmentation 

nucleation  and  condensation 
vaporization,  sublimation,  and  melting 

Replacement  of  TBL  -  integral  or  finite  difference  method? 

mass  diffusion? 
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DESIGN  COMPARISON 


STAR  37X  ROCKET  MOTOR 


Y 107 9099  /  elkton  division 


Y1079100  //t^O^COC  /  ELKTON  DIVISION 


APPARATUS  FOR  CROSS  FLOW  PHOTOGRAPHY 


HIGH-SPEED 

CAMERA 


CROSS-SECTION  THROUGH  WINDOW 
ROTOR  SHOWING  NORMAL  TO  THE 
BURNING  SURFACE  CONFIGURATION 


A)  ARRANGEMENT  FOR  PHOTOGRAPHING  ALUMINIZED  PROPELLANT. 


FLOW  IN  PORT 
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C)  PROPELLANT  CONFIGURATION 
FOR  NORMAL  TO  THE  BURNING 
SURFACE  PH0T06RAPHT. 


B)  PROPELLANT  CONFIGURATION 
OF  SIDE  VIEW  PHOTOGRAPHY. 


PARTICLE  DISPLACEMENT  VS  TIME 


TIME,SECONDS(*/0?> 
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Appendix  15:  Expended  Inerts  and  Slag  Performance  Procedure 
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J.  T.  LAMBERTY 
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e)  Time  interval  used  for  impulse  determination 

f)  Propellant  mass  assumption 
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Appendix  17:  Particle  Impingement  Modelling  -  Performance  Effects 


I 

I 

I 

l 

1 

1 

I 

1 

I 

1 

I 

1 


229 


STATE-OF-THE-ART  MOTORS 


MOTOR  EFFICIENCY 
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Appendix  20:  Efficiency  Definitions 
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Appendix  22:  Unsteady  Three-Dimensional  Flow  in  Propulsive  Nozzles 


UNSTEADY  THREE-DIMENSIONAL  FLOW 
IN  PROPULSIVE  N0Z2LES 


Joe  D.  Hoffman 
Purdue  University 


February  1980 
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Combined  Nozzle  and  Missile  Rotation,  Figure  2,  continued 


Gas  Dynamic  Model 


1.  Isentropic  (i.e.,  invtscid,  adiabatic)  flow. 

2.  No  body  forces  except  those  due  to  missile  rotation. 

3.  No  condensed  phases. 

4.  Perfect  gas. 

5.  Frozen  flow. 


817 


1.  Unsteady  three-dimensional  flow. 

2.  Steady  flow  as  limit  at  large  time. 

3.  Cylindrical  coordinate  system. 

4.  Axi symmetric  boundaries. 

5.  Transformed  computational  space. 

6.  Three-dimensional  effects  due  to  rotation  only. 

7.  MacCormack's  method  at  interior  points. 

8.  Kentzer's  method  at  all  boundary  points. 
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(a)  Rocket  motor  flow  geometry. 


Solid  Boundary 


Boundary  Points 

(b)  Propulsive  nozzle  unit  processes. 


Figure  5.  Flow  field  definition. 


Bivariate  Interpolation  Employed 
to  Obtain  the  Cylindrical  Properties 


Linear  Interpolation  Employed  to 
Obtain  Cartesian  Properties  at 
Points  9, 10,  II,  and  12. 


Figure  9.  Cartesian  grid  network  in  the 
region  of  the  centerline. 


A**.-*"*  ■ 


Propulsive  nozzle  rotation  study  geometry 


TABLE  1.  Data  for  the  nozzle  studies. 


Flow  Properties 


Stagnation  Pressure,  psia 

4570 

Stagnation  Temperature,  R 

5987 

Ratio  of  Specific  Heats 

1.15 

Molecular  Weight,  products 

26.71 

Gas  Constant,  f t-1 bf/ ( 1 bm-R) 

57.85 

Ambient  Pressure,  psia 

14.7 

Propellant  Data 

Burn  Rate  at  2000  psi,  in. /sec  11.0 

Pressure  Exponent  0.7 

3 

Density,  1  bin/ in 


0.063 


(c)  Velocity  field  in  the  region  of  the  centerline 

800  deg/sec. 

Figure  16,  continued. 


Figure  17.  Exit  plane  pressure  distortion. 
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(b)  Missile  rotation,  400  deg  /sec. 


400  dps 


Figure  17,  continued. 


200  dps  iso 

(c)  Missile  rotation, 


*1^*400  dps 
*800  deg/sec. 


Figure  17,  continued. 


Torque.in-lbfxiO"* 

-34-  _ 

Missile  Rotation* 800  deg/sec 


1  + 


Missile  Rotation  *400  deg/sec 


Nonrotating  Missile 


200  400 


Nozzle  Rotation 
deg/sec 


2 


3t 

(a)  Torque. 


Effects  of  rotation  on  the  forces  and 
moments  acting  on  the  nozzle. 


Side  Force,  Ibf  xlO“ 


Figure  14.  continued. 


400  -  200 


200 

400  Rotation 

deg /sec 

Symbol 

Missile  Rotation 

o 

0 

o 

400  deg /sec 

A 

800  deg/sec 

(c)  Axial  force. 


Figure  14,  continued. 
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Appendix  23:  Advanced  Performance  Prediction  Methodology 


ADVANCED  PERFORMANCE  PREDICTION  METHODOLOGY 
JAY  N.  LEVINE 

The  goal  of  this  discussion  period  will  be  to  identify: 

1.  Limitations  of  present  performance  prediction  methodologies 

a.  SPP 

b.  Technology  In  general 

2.  Impact  of  solid  rocket  motor  technology  developments  on  performance 
prediction  methods. 

The  following  items  are  listed  to  stimulate  thinking  about  the  subject. 
These  items,  as  well  as  others  suggested  at  the  meeting,  will  be  discussed. 

An  effort  will  be  made  to  catergori?e  the  items  into  three  groups. 

1.  Those  items  which  are,  or  are  very  likely  to  be.  important  enought  to 
definitely  warrant  that  the  ability  to  treat  them  be  included  in  an  advanced 
performance  prediction  methodology. 

2.  Those  items  which  might  be  Important,  and  should  be  carefully  con¬ 
sidered  for  Inclusion  in  an  advanced  performance  prediction  methodology. 

3.  Those  items  which  are  not  likely  to  be  important,  or  have  a  significant 
effect  on  performance. 
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1.  Effect  of  all  types  of  thrust  vector  control  methods  on  performance 


a.  Canted  nozzles 

b.  Glmbaled  nozzles 

c.  Jet  interaction  TVC 

d.  Boundary  layer  TVC 

2.  Grain  asymmetry,  multidimensional  ballistics 

3.  Non-circular  cross-section  nozzles 

4.  Scarfed  nozzles 

5.  ENEC  and  very  high  area  ratio  nozzles 

6.  Spin  effects 

7.  Nozzle  wall  roughness  and  asymmetrical  nozzle  erosion 
•8.  Exotic  propellant  ingredients 

9.  Multi -propellant  motors 

10.  Aerodynamic  heating  effects 

11.  Extreme  motor  operating  conditions 

In  regards  to  the  various  loss  mechanisms: 

Two-phase 

Divergence 

Boundary  layer 

Kinetics 

Submergence 

Erosion 

Combustion 

Mass  transfer  between  phases 
Impingement 

Which  do  you  feel  are  modeled  satisfactorily,  at  present?  Which  do  you  feel 
require  additional  work? 
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Results 

The  old  familiar  and  more  recent  drag  models  have  been  compared 
to  the  data  of  UTC  (Ref.  1),  Zarin  (Ref.  8)  and  Bailey/^liatt  (Ref.  9 
for  Rep  £  40) .  For  each  of  the  51  UTC  data  points ,  46  Zarin  data 
points,  and  45  Bailey  /Hiatt  data  points  the  predicted  Cp  was  normal¬ 
ized  by  the  measured  Cn.  An  average  and  root-mean-square  (rms)  error 
was  then  calculated  for  the  j—  set  of  data  from 


A  small  rms  error  implies  that  the  model  correlates  the  data  with 
very  little  scatter,  while  a  small  average  error  implies  that  the 
model  predicts  values  near  the  middle  of  the  scatter-band.  The 


results  are  presented  in  Table  III  and  the  rms  errors  are  summarized 


below : 


UTC 


Nitrogen 

Freon 

(35  points) 

(l6  points) 

Kliegel  (SPP) 

2.17 

1.65 

Crowe  (TPP) 

0.14 

0.16 

Walsh 

0.49 

0.33 

Henderson 

0.15 

0.20 

Hermsen 

0.11 

0.11 

Zarin 

Bailey/Hiatt 

All 

y  £  1.77 

p 

U  >  1 .77 
P 

(46  points) 

(12  points) 

(32  points; 

0.09 

0.09 

0.24 

0.21 

0.27 

0.15 

0.06 

0.02 

0.26 

0.09 

0.05 

0.05 

0.06 

0.10 

0.30 

The  data  is  plotted  together  in  Fig .  4 . 


Note  that  the  Bailey /ttiatt  data  was  arbitrarily  divided  into  two  groups 
by  Mp  —  1.77  (i.e.  at  the  maximum  Mp  in  the  UTC  data  at  which  the 
Hermsen  model  was  correlated).  It  is  unlikely  that  Mp  >  1.77  will 
ever  be  encountered  in  a  rocket  nozzle  (although  it  might  be  In  plumes) 
because 


1)  ■*  4(r  shows  that  Id  < 

p  g  p 

region  of  the  nozzle,  even 

particle  (aV  -  ug) . 


1.0  everywhere  in  the  subsonic 
for  the  worst  case  of  a  stationary 
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2)  Mp *** yspj  0.7  for  a  mean-size  particle  in  a  53:1  nozzle 

(see  Fig.  1)  . 

The  significance  of  the  error  in  predicting  the  drag  coefficients 
at  high  Up  is  limited  by  the  fact  that  there  was  a  data  scatter  of 

nearly  20%  in  high-U  low-Re  test  results. 

P  P 

Finally,  all  the  above  data  was  obtained  in  the  absence  of  gas- 
stream  turbulence.  Some  turbulence  should  exist  in  rocket  motors, 
but  Zarin  has  shown  that  the  effect  is  probably  negligible  at  the  low- 

Re  encountered  there . 

P 

Conclusions 

The  above  table  shows  that: 

1)  As  expected  the  Kliegel  model  does  well  in  the  low-Mp 

si ip /transit ion  regime  represented  by  the  Zarin  and  Bailey/ 
Hiatt  data  but  fails  badly  as  the  free-molecular  regime  is 
encountered  in  the  UTC  data . 

2)  As  expected  the  Crowe  empirical  relation  does  well  for  the 
free-molecular  conditions  of  the  UTC  data  it  was  forced  to 
fit,  but  becomes  less  accurate  for  the  slip/transition  regime 
of  the  Zarin  and  Bailey  /Hiatt  data  that  it  was  not  forced 

to  fit. 

3)  As  expected  the  Walsh  correlation  does  well  in  fitting 
the  Zarin  and  low-Mp  Bailey /Hiatt  data  as  it  was  forced 
to  do,  but  fails  to  predict  the  free-molecular  UTC  data 
or  the  high-Mp  Bailey /Hiatt  data. 

4)  The  Henderson  model  does  quite  well  in  predicting  all  the 
data  at  -  1.40,  but  not  so  well  for  freon  (  Y  =*  1.09). 

5)  The  Hermsen  model  is  the  best  predictive  technique  for 

Mp  is  1.77  with  small  average  and  rms  error  for  both  free- 
molecular  and  slip-flow  regimes;  note  that  the  Hermsen  rms 
error  in  predicting  the  UTC  data  is  no  worse  that  the  re¬ 
ported  experimental  accuracy  (average  of  +  11%  for  the  51 
points  -  see  Table  I) . 

6)  Thus  the  Hermsen  model  should  be  used  for  nozzle  flows 
while  the  Henderson  model  may  be  better  for  plume  flows 
where  Mp  might  exceed  1.77  and  Yg  is  approaching  1.40. 
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Appendix  25:  Problem  -  Measurement  of  Thrust  for  a  Solid  Rocket 
Fired  in  the  Vertical  Attitude 
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Major  sources  of  error: 
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(301)  953-7100,  x3151 _ 


(805)  277-5127 


NASA,  Marshall  Space  Flight  Center 
Attn:  Mr.  B.  W.  Shackelford 
Code  EP25 

Huntsville,  AL  35812 
(205)  453-3808 _ 


AFRPL  (MKPA) 

Attn:  Dr.  F.  Roberto 
Edwards  AFB,  CA  93523 


ARO,  Incorporated  (AEDC) 

Attn:  R.  Bruce  Runyan,  ETF-J 
Tullahoma ,  TN  37389 

(615)  455-2611,  x7373 _ 

Thiokol  Corporation 
Wasatch  Division 
Attn:  Dr.  M.  Salita , IAS  ZRl 
PO  Box  524 

Brigham  City,  UT  84302 
(801)  863-2163 _ 


McDonnell-Douglas  Astro.  Co/East 
Attn,:  Mr.  Joseph  Smuckler,  Dept.  E243 
PO  Box  516 
St.  Louis,  MO  63166 

(314)  232-7577 _ 


Rockwell  International /Rocketdyne 
Attn:  Mr.  George  Sopp,  D/545-113,  AC06 
6633  Canoga  Avenue 
Canoga  Park,  CA  91304 

(213)  884-3446 


’ASA,  Headquarters 

Attn:  Mr.  F.  W.  Stephenson,  Jr.,  RP 
600  Independence  Avenue  SW,  Room  625 
Washington,  DC  20546 

(202)  755-3280 _ 


jet  Liquid  Rocket  Company 
»"n  *»t.  J.  W.  Salmon 

a...  M77? 

».  .  .  CA  95813 


Naval  Weapons  Center 
Attn:  Dr.  Charles  J.  Thelen 
Code  3205 

China  Lake,  CA  93555 
(714)  939-7206 _ 
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NASA,  Johnson  Space  Center 
Attn:  Mr.  J.  G.  Thibodaux,  BP 
Houston,  TX  77058 


ARO,  Incorporated  (AEDC) 

Attn:  Floyd  E.  Turner,  Jr.  ETF-J-EB8 
Tullahoma,  TN  37389 

(615)  455-2611,  x452 _ 


TRW  Systems  Group 
Attn:  Mr.  T.  J.  Walsh 
Bldg.  523,  Dept.  315 
PO  Box  1310 

San  Bernardino,  CA  92402 
(714)  382-6089 _ 


Rockwell  International  Corporation 

Rocketdyne  Division 

Attn:  Mr.  W.  T.  Webber 

6633  Canoga  Avenue 

Canoga  Park,  CA  91304 


AFRPL  (CA) 

Attn:  Dr.  R.  R.  Weiss 
Edwards  AFB,  CA  93523 


(805)  277-5622 


Thiokol  Corporation 

Huntsville  Division 

Attn:  Mr.  R.  H.  Whitesides ,7r  ,81^. 

Redstone  Arsenal 

Huntsville,  AL  35807 

(205)  876-3628 _ 


t 
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PUB  312  DISTI IBOTXOB 


-  ARBI  - 

Aray  Bsl  coad/Redstone  Arsenal 
1  Xfahin,  P.  J. 

1  Radke,  R.  P. 

-  RAVI  - 

Ratal  ord  sta/Indian  Haad 
1  Johnson,  C, 

1  Lyon,  H. 

Ratal  Hpn  Ctr/China  Lake 
1  Derr,  8.  L. 

1  Eacalllar,  P.  n. 

1  Pabans,  8, 

1  Tbelao,  c#  J« 

-  AIR  roRCE  - 

A  Fi<  PL/Ed  wards  APB 
1  Beckaan,  C. 

1  Gelsler,  R.  L. 

1  George,  D. 

1  Latin*,  J. 

1  Roberto,  P. 

1  Selph,  C.  C. 

1  Weiss,  B.  R. 

-  RASA  - 

RASA  -JSC/Houston 
1  Kahl,  R.  C. 

1  Thibodaax,  J.  G. 

RASA -La  RC/Haapton 
1  Buahnell,  D. 

RASA -LeRC/Cl eta land 
1  Blttker,  D. 

1  Gordon,  3. 

NASA-Bgat  Ops/Washington,  DC 
1  Stephenson,  P.  a. 

RASA -RSPC/H arshal 1  Space  Plight 
Center 

1  Forsythe,  D.  J. 

1  Gross.  K. 

1  Richaond,  R. 

1  Shackelford,  B.  R. 

-  liOR-GOVERNBERT  - 

Acurer/Roantain  flee 
1  Abbott,  B. 

1  Chiba,  Z. 

1  Daha,  T. 

1  Etana,  R.  B. 

A-rojet  Liquid  Propulsion/ 
Sacraaento 
1  Ito,  J,  X • 

1  Pieper,  J,  L. 

1  Salaon,  J,  a. 

Aerojet  Solid  Propulsion/ 
Sacraaento 
1  Chan,  G.  0. 

1  Ditore,  B.  J, 

1  Bockeahaupt,  J.  o. 

1  Seibert,  j.  R. 
Aeroapace/Los  Angales 
1  Chang,  I-Shih 

1  Landsbaan,  E.  B, 

ARO/Tullaboae 

1  Little.  R.  R. 

1  innyan,  I,  B. 

1  Turner,  Jr.  ,  p,  ?. 

Atlantic  isch/Aleiandria 
1  Headaraon,  r, 

1  King,  B. 

1  Procinaky,  I,  n. 

Roe  a  ng/Saattle 

1  Groan,  a.  I, 


-  R  OR -GO  VP.  PH  BERT  -  (con**d) 

(con t'd) 

Boeing/Seattle  (cont*d) 

1  Kernes,  J. 

1  BcAndrews,  E. 

California  Inst  lee h-JPL/ 
Pasadena 

1  Anderson,  P.  A. 

Cohen  Prof  Stcs/hedlands 
1  Cohen,  R. 
Continuua/Huntstille 
1  Prozan,  R.  J. 

CPIA/ Laurel 

1  Hege,  H.  F. 

1  Reedy,  T.  L. 

Heccules/R  agna 

1  Boyd,  D.  L. 

1  Daines,  H.  L. 
Hercules/McGregor 
1  Brooks,  R .  T. 

KVB/Tustin 

1  cherry,  S. 

1  Hunter,  S.  c. 

1  Kliegel,  J.  R. 

Lockheed  Bsl  S  Space/Huntstille 
1  Preeaan,  J.  A. 

1  Penny,  B. 

1  Snitb,  S.  D. 

Lockheed  Bsl  6  Space/Sunny tale 
1  Hoosaan,  R. 

LTV  Aerospace/Baapton 
1  Bahn,  G.  S. 

Barquardt/Van  Buys 

1  Caapbell,  J.  G. 

BcDonnell  Douglas/St.  Louis 
1  Sauckler,  J. 

Purdue  Onit/Weat  Lafayette 

1  Hoff aan,  J. 
Rockeall/Canoga  Park 
1  Ganiac,  V. 

1  Sopp,  G. 

1  Habber,  H.  T. 

Scianca  Applications/St.  Louis 
1  Gold  ford,  k. 

Software  S  Eng  Assocs/Santa  Ana 
1  Coats,  D.  E. 

1  Eickerson,  G»  R. 
Thiokol/Br ighaa  City 
1  kbel,  R. 

1  Lesley,  G. 

1  Salita,  B. 

1  Sayar,  L. 

Thiokol/Elkton 

1  Brundige,  u.  N. 

1  Kirnchner,  Z. 
Thlokol/Huntatllle 

1  Blkkelsen,  C.  D. 

1  Ovens,  T.  F. 

1  Whitesides,  Jr.  ,  B.  H. 
TRW/Redondo  Beach 
1  Harry,  D. 

TRV/San  Bernardino 
1  Ralsh,  T.  J. 

United  Technologies-CSD/Sunny- 
tale 

1  Harasan,  8,  H. 

1  Laaberty,  J.  T. 

United  Tachnologies-Ptw/west 
Pals  Beach 
1  Carroll,  R.  G. 

1  Bayes,  T.  C. 

Unit  Utah/Sait  Lake  City 
1  flandro,  G,  A. 

Washington  State  Unit/Pullaan 
1  Crowe,  C.  T. 
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